Athermodynamic Alignment in the Two-Dimensional Enstrophy Cascade 
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We study inertial-range statistics in the direct enstrophy cascade of two-dimensional turbulence 
via a numerical simulation of the forced Navier-Stokes equation. In particular, we obtain the dis- 
tribution of the enstrophy flux and of the angle of alignment of the large-scale vorticity gradient 
with the small-scale vorticity transport vector. These distributions are surprisingly symmetrical 
and "athermodynamic" , not explainable by a local eddy- viscosity approximation with coefflcient of 
either positive or negative sign. By a systematic evaluation of the role of various triadic interactions, 
we can trace the origin of the strong symmetry to the dominance of infrared non-local triads in the 
enstrophy cascade. 
PACS numbers: 47.27.Ak,47.27.Gs 



Two-dimensional (2D) turbulence has been a fascinat- 
ing topic for over thirty years, since the seminal papers of 
Kraichnan [Q , Batchelor |^ and Leith Q . The enstrophy 
cascade to small-scales proposed in those works is the 2D 
analogue of the energy cascade in three-dimensional (3D) 
turbulence. This forward cascade of enstrophy continues 
to be the subject of many recent theoretical and 
experimental investigations. One of the reasons for 
the sustained interest in the subject is its importance for 
the interpretation and analysis of atmospheric dynamics. 
In addition, it is topic of great fundamental interest in 
statistical physics. For example, in free turbulent decay 
a flux to high wavenumber spontaneously develops for 
the inviscidly conserved variables, energy in 3D and en- 
strophy in 2D. This is a striking example of irreversible 
behavior arising from a time-reversible dynamics, the 2D 
and 3D Euler equations. Recent research on 3D tur- 
bulence has identified the key geometric and statistical 
property behind this forward cascade as certain inertial- 
range alignments between characteristic large-scale and 
small-scale quantities In the opinion of the au- 

thors, one of the most crucial goals of modern turbulence 
theory is to understand how such alignments are spon- 
taneously produced by essentially inviscid dynamics at 
high Reynolds number. In this letter we study the cor- 
responding alignments in 2D, which lead to the forward 
enstrophy cascade. The main new result of our study 
is that the 2D alignments are surprisingly weak, much 
more so than in the 3D energy cascade. In a sense that 
shall be explained below, the most likely alignment in 
2D is athermodynamic, leading neither to forward nor to 
backward cascade. Any successful theory of 2D turbu- 
lence must account for this remarkably time-symmetric 
behavior. We shall establish a relation with another well- 
known property of 2D enstrophy cascade — its domina- 
tion by nonlocal triadic interactions |^,^ — which distin- 
guishes it sharply from the local energy cascade in 3D. 
Our results have important consequences for the numer- 
ical modeling of large-scales in turbulent 2D flows. 

We have simulated the equation 



dtiJ + vVw + Vi{-^)-P'Lo + z/„(-A)P"cj = F (1) 

in a square domain with side L — 27: and periodic bound- 
ary conditions. Here v is the velocity and a; = V x v is 
the vorticity, is a stirring force applied to wavenumbers 
k| = 4 — 7 to give a constant enstrophy input rate. We 
add hyperviscosity with p„ = 8 at high wave numbers 
to extend the inertial range. We also add hypoviscosity 
with Pi = 2 at small wave numbers to destroy box-size 
vortices. The equation was solved using a pseudospectral 
parallel code with full dealiasing. The time stepping was 
a second-order Adams-Bashforth method. The resolution 
was 2048^. A statistical stationary state was achieved af- 
ter evolved for about 200 large-eddy turn-over time. In 





Fig.l we plot a visualization of the vorticity field of our 
flow at one instant. It separates into two distinct re- 
gions: large-scale, coherent vortex structures and strain- 
dominated regions between these, filled with drawn-out, 
fine filaments of vorticity. It is in the latter regions where 
the enstrophy cascade occurs. In Fig. 2 we plot the spec- 
tral enstrophy flux as a function of wavenumber. Clearly, 
we have about a decade and a half of inertial range, where 
this flux is constant. The inset of Fig. 2 shows the energy 
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spectrum of this final steady-state. The power-law range 
is slightly steeper than the —3 law predicted in [|UE|, but 
consistent with the later logarithmic correction |12[| . 
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FIG. 2. Enstrophy flux spectrum, rj is enstrophy dissipa- 
tion. Inside box is energy spectrum. 

We wish to study the statistics of the enstrophy cas- 
cade locally in space. We therefore consider a local flux, 
which measures transfer of enstrophy into small scales 
at a fixed point in physical space. The theoretical inter- 
est of such a quantity for studying the intermittency of 
the 3D energy cascade was first emphasized by Kraichnan 
1^ , who used banded Fourier series in his definition. We 
use instead a smooth filter to differentiate the large-scale 
and small-scale modes. This is the same method used in 
the large-eddy simulation (LES) modelling scheme p^ ] 
and in our earlier discussion of the 3D case We 
apply the filter to the 2D Euler equations in their vortic- 
ity formulation dtOJ -I- (v • V)w — 0. That is, we consider 
the "large-scale vorticity" defined as the convolution field 
uji = Gi * oj, with some suitable filter function G^, and 
large-scale velocity field likewise defined by V£ = * v. 
In particular, in this paper G^ is the Fourier space Gaus- 
sian filter. The equation obtained by low-pass filtering 
is 

dt-UJi{Y,t) + V •[wt{Y,t)Wi{r,t) + (Ti{r,t)] ^ Q. (2) 



Here cri = (vw)^— v^ZU^ is a vector representing the space- 
transport of vorticity due to the eliminated small-scale 
turbulence. It plays the same role in 2D as the turbu- 
lent stress tensor ti in the analogous 3D equation for 
the large-scale velocity. From the previous equation, a 
balance equation is easily derived for the local density 
he{r,t) = ilU^(r,t) of the large-scale enstrophy: 



dthi{r, t) + V-K^(r, t) = -Zf(r, t) 



(3) 



in which the current Kf(r, i) = h£{r,t)'ve{r,t) + 
uJi{r,t)cre{r,t) represents space-transport of the large- 
scale enstrophy, and 



Ze{r,t) = -VLJi{r,t)-(7e{r,t) 



(4) 



irreversible forward cascade of enstrophy occurs precisely 
when it develops a positive mean-value. Notice that this 
positivity is equivalent to the thermodynamically natu- 
ral statement that the turbulent vorticity transport (Ti 
should tend to be "down-gradient" , that is, anti-parallel 
to the large-scale vorticity gradient Vw^. The required 
statistical anti-correlation of cr^ and VoJ^ is an alignment 
property characteristic of the 2D enstrophy cascade. It is 
analogous to the much-studied alignment of tensor quan- 
tities in 3D, the stress-tensor due to small-scales and 
the large-scale strain S^, which underlies the energy cas- 
cade to high-wavenumbers | |ll| . 

One of the major results of our numerical study is the 
probability density function (PDF) of the enstrophy flux 
in the steady-state cascade, P{Zi), shown in Fig. 4 for 
several filtering length-scales £ in the forward cascade 
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FIG. 3. Normalized PDF of enstrophy flux Ze{r,t) at dif- 
ferent fllter length £. 




is the enstrophy flux out of the large-scales into the small- 
scale modes. This quantity is odd under time-reversal; an 
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FIG. 4. PDF(6'), 9 is the angle between ere and VZU^ at 
£ = 7r/130. 

range. The analogous PDF of energy flux has been ob- 
tained numerically in |l6| , but, to our knowledge, this is 
the first result for enstrophy flux in 2D. The most strik- 
ing fact about this PDF is its near symmetry. In nearly 
50% of the realizations, the enstrophy flux is negative, or 
backward to large-scales. The skewness of the PDF is also 
quite small, only 0.15. In contrast, the PDF of energy 
flux in 3D has only about 33% of the realizations with 
negative values and a skewness of 11 in the inertial range 
[ pTf . The characteristic feature of the energy flux PDF in 
3D is the asymmetrically long tail to the right, whereas 
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for the enstrophy flux in 2D we see that tails to the right 
and left are about equal. Surprisingly, enstrophy in 2D 
is quite reluctant to cascade to high wavenumber. To 
further quantify this, we show in Fig. 4 the PDF of the 
angle of alignment 9 between the vectors cr^ and Va;^. 
We consider only the filtering length £ — vr/lSO in the 
inertial-range, because the results vary little with £ in 
that range. Consistent with our result on the PDF of 
the flux, we see that the distribution is nearly flat, which 
in 2D is characteristic of a pair of independent, isotrop- 
ically distributed vectors. There is a slight tendency for 
the two vectors to be anti-correlated but it is very small: 
their correlation coefficient is only = —2.53 x 10"'^! In 
fact, we see from Fig. 4 that the (weakly) most probable 
configuration of the pair of vectors is to be perpendic- 
ular. That is, there is more of a tendency of vorticity 
transport cri to be along streamlines of the large-scale 
fiow than to be down the vorticity gradient. Such an or- 
thogonal alignment of the vectors is "athermodynamic" , 
producing neither forward cascade as does the thermody- 
namic, down-gradient alignment, nor backward cascade 
as does an anti-thermodynamic, up-gradient alignment. 
The entire forward cascade phenomenon in 2D is only a 
weak, secondary effect. 

To shed further light on this, we consider a sequence 
of successive approximations to cr^, in order to study the 
effect of neglected interactions on the enstrophy cascade. 
The first approximation we make is the 2D analogue of 
the Similarity Model in 3D . The basic assumption of 
this model is that the wavevector triads with (at least) 
one mode at scales ^ I do not contribute much to cr^. 
To implement this idea, let us define a new lengthscale 
t < I and low-pass filter v of v with scales < £ removed. 
Then, the Similarity Model is 

rrf' ^CsM^-^^ (5) 

where Csm is a constant. In fact, it is reasonable to take 
CsM ~ 1 and £ = £, which hereafter we do. The as- 
sumption of the Similarity Model is generally held to be 
well-satisfied, in 2D as well as in 3D. The reason is that 
the scales < i are believed to make only a disorganized, 
uncoordinated contribution to the vorticity transport cri. 
This means that contributions of those scales are subject 
to large cancellations, producing an extra small factor 
£. See [|2|. Thus, we believe that the Similarity Model 
omits only contributions that are already negligible in 
reality. This expectation is borne out by our simulation. 
In fact, we have found that cr^ and cr^^ are correlated at 
the 99% level, which is a higher level of correlation than 
the model enjoys even in 3D. The PDF of the enstrophy 
fiux in the Similarity Model is so similar to that for the 
true enstrophy fiux in Fig. 3 that it will not be shown 
here; in particular, it retains the high symmetry of the 
exact flux. Although these are only "a priori" tests of 
the model, we expect that the Similarity Model is likely 
to perform well also "a posteriori" , in an actual 2D LES 
computation, at least in a "mixed" version with addition 
of a small eddy- viscosity. 



Although the Similarity Model neglects distant small- 
scales, it retains infrared (IR) nonlocal triadic interac- 
tions. The latter are believed to be the crucial inter- 
actions in the 2D enstrophy cascade 0,^, which pro- 
ceeds by stretching of small-scale vorticity-gradients due 
to strain arising from the largest-scale vortices. These IR 
nonlocal triadic interactions lead to an enstrophy transfer 
which is ultra-local in wavenumber: transport is diffusive 
in fc-space, a (biased) random walk by very small steps 
with size proportional to the small wavenumber of the 
largest-scale vortices In particular, Kraichnan 

has argued |l8) that such ultra-locality of the enstrophy 
cascade in wavenumber space will prevent the validity of 
any model for erg which is local in physical space. In 
fact, the expression erf^'^ in Eq.(||) satisfies Kraichnan's 
criterion: it is spatially non-local, given by an integral 
(convolution) over length-scales of order £. 

To test the importance of these IR nonlocal triadic 
interactions, we can make a local expansion of the Sim- 
ilarity Model in physical space. The leading-order term 
is the analogue of the so-called Nonlinear Model in 3D 
0: 

erf ^ = C2£^T>rV[;}i, (6) 
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FIG. 5. PDF(6l), 6 is the angle between r^"^ and Vue at 
different filtering length £. 

where Dij = dvi/dxj is the large-scale velocity gradient 
(or deformation) tensor and C2 is the 2nd-moment of the 
filtering function G. Neglected higher-order terms are 
formally smaller by additional powers of £. This expres- 
sion allows for a very intuitive understanding of the ten- 
dency for enstrophy to cascade forward to small scales [|j . 
In fact, the compression of vorticity level sets close to- 
gether in the strain-dominated regions outside the large- 
scale vortices should tend to produce strong gradients of 
the vorticity along the compressive direction. Such an 
alignment has long been known to occur in the 2D dissi- 
pation range, between the full vorticity gradient Vlj and 
r^"', the right eigenvector of D with negative eigenvalue 
|]l9| . See also j2^. We have verified that the same re- 
sult holds for the inertial-range quantities, VcD and D^, 
defined by filtering. In Fig. 5 is shown the PDF of the 
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angle between and Vw£ for various filter lengths t 
in the inertial-range. We see a strong tendency of these 
two vectors to be parallel, especially at the smaller scales. 
In consequence, we can roughly approximate cr^^ in the 
strain regions as 

(xf ^ « -C2eb,VZ:,, (7) 

where —D^ is the negative eigenvalue of the deformation 
matrix at length-scale I. Note that Dg ^ D for all £, 
because it is dominated by the largest scales. We have 
obtained finally an eddy-viscosity approximation, with 
turbulent viscosity vx C2DP , very similar to that first 
proposed by Leith ||^ . See also . 

However, the considerations advanced by Kraichnan 
10 have cast doubt a priori on the validity of any such 
approximation local in physical space. We have seen 
earlier that the prediction of down-gradient transport 
yielded by Eq.(|^) is very poorly satisfied in our simu- 
lation. In fact, the full Nonlinear Model in Eq.(^) sub- 
stantially overpredicts the forward cascade of enstrophy. 
In Fig. 6 we plot the PDF of the Nonlinear Model of en- 
strophy flux, Z^^ = -C2i'^{Vu))Jl)(Vu}i, calculated in 
our simulation. We see that the model prediction is quite 
similar to the energy flux PDF in 3D, with a strong skew- 
ness '^8.1 and a long, asymmetrical tail to the right. 
However, we have already found that such a prediction 
for the 2D enstrophy cascade is totally spurious and that 
the true enstrophy flux is distributed nearly symmetri- 
cally over positive and negative values. The correlation 
coefficient of cr^^ with the true cr^ is only 0.03 in our 
simulation, for £ = 7r/130. This value may be contrasted 
with the 3D case, where the correlation coefficients be- 
tween components of the true stress tensor t^j, and the 
nonlinear model r/^^ are of the order of 90% pl . 
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FIG. 6. Normalized PDF of Nonlinear Model of enstrophy 
flux at ^ = 7r/300. 

We conclude that the statistics of enstrophy flux in 2D 
are far more symmetrical than those of energy flux in 



3D. Correspondingly, the alignment between the large- 
scale "force" VaJ^ and the small-scale "flux" cr^ is much 
weaker than is the alignment of corresponding quantities 
in 3D. To the extent that there is any alignment at all, it 
is mainly "athermodynamic" . Any final theory of 2D tur- 
bulence must account for this striking difference from the 
3D forward cascade. While offering no complete explana- 
tion, we have shown that the phenomenon is connected 
to the dominance of nonlocal triadic interactions in 2D, 
whereas energy cascade in 3D is predominantly by local 
triads. In consequence, a spatially-local approximation 
to the 2D enstrophy flux fails, whereas the correspond- 
ing approximation in 3D correlates very well with the 
facts. 
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